/FILN,PRESSURE_VESSEL
/TITLE,THE OPTIMAL DESIGN OF TANK   
!ͼһµķű
*SET,R1,1508                   ! Բڰ뾶        
*SET,R2,1558                   ! ڰ뾶     
*SET,T1,200                    ! ԲǱں
*SET,T2,100                    ! Ǳں
*SET,H,1600                    ! ߶
*SET,L,300                     ! ׶ζγ
*SET,PRESS,10                  ! ѹ
*SET,THEROY,PRESS*R2/(2*T2)    ! ͲӦֵ

/PREP7  
ET,1,PLANE82 
KEYOPT,1,3,1                   !ΪԳ
MP,EX,1,1.765E5                ! ģ
MP,PRXY,1,0.3                  !ɱ
*AFUN,DEG                      !ĽǶȵλ
ANGLE=ASIN(L/R2)               !׶ζǵ
!ģ
BLC4,R1,-(H-L),T1,(H-L)         
PCIRC,(R2+T2),R2,90,ANGLE
A,4,3,5,8
AADD,ALL

!񻮷
ESIZE,35,0                     !ûֵԪС 
MSHAPE,0,2D                    !ıл
MSHKEY,0                       !趨ɷ
AMESH,ALL                      !

!ýڵţΪ·
CSYS,1
NSEL,S,LOC,X,R2                
NSEL,R,LOC,Y,ANGLE
*GET,NUM_1,NODE,0,NUM,MAX

!׶ζת۴ؾھNODEţλΪ(R2)
NSEL,S,LOC,X,R2+T2
NSEL,R,LOC,Y,ANGLE
*GET,NUM_2,NODE,0,NUM,MAX   
!׶ζת۴ؾ⾶NODEţλΪ(R2+T2)
FINISH  

/SOLU   
CSYS,0
NSEL,S,LOC,Y,-(H-L)
D,ALL,UY               !ģ¶ʩYԼ
CSYS,0
NSEL,S,LOC,X,0
DSYM,SYMM,X            ! ģX=0㴦ڵʩXԳԼ
LSEL,S,LINE,,4,10,3
NSLL,S,1
SF,ALL,PRES,PRESS                     !ڱڽڵʩӾѹ
ALLSEL  
SAVE
SOLVE                                 !
SAVE,PRESSURE_VESSEL_RESULT           !DBļ
FINISH  

/POST1  
PATH,PATH1,2,30,20                     !·
PPATH,1,NUM_1                          
PPATH,2,NUM_2
PDEF,MISES,S,EQV,AVG                   ! ӳӦ·
*GET,MISES,PATH,,MAX,MISES             ! õ·Ӧֵ
INTENS2=MISES/THEROY                   ! Ӧϵ
LGWRITE,PRESSURE_VESSEL,LGW,,COMMENT   ! Żļ
FINISH  

/OPT                                   ! Ż
OPANL,'PRESSURE_VESSEL','LGW'          ! ȡŻļ
OPVAR,L,DV,290,320                     ! Ʊ
OPVAR,ANGLE,SV,0,30                    ! ָ״̬
OPLOOP,PREP,PROC,SCALAR
OPVAR,INTENS2,OBJ,,,0.0001             ! ָĿ꺯INTENS2
OPSAVE,'PRESSURE_VESSEL','OPT'         ! Ż״̬
SAVE,PRESSURE_VESSEL_MODEL    
OPTYPE,FIRS                            ! ָŻΪһ׷
OPFRST,20,,,                           ! ָѭ
OPEXE                                  ! ִŻ
OPKEEP,ON                              ! 
OPLIST,ALL,,0                          ! бʾѡ
SAVE,PRESSURE_VESSEL_OPT_RESULT           
/AXLAB,X,ITERATION NUMBER              ! 趨X
/AXLAB,Y,STRUCTURAL INTENS             ! 趨Y
PLVAROPT,INTENS2                       ! Ŀ꺯ѭ仯ʾ
/AXLAB,X,L                             ! 趨X
/AXLAB,Y,INTENS                        ! 趨Y
PLVAROPT,L,INTENS2                     ! Ŀ꺯L仯ʾ
FINISH


 LIST OPTIMIZATION SETS FROM SET   1 TO SET  21 AND SHOW
 ONLY OPTIMIZATION PARAMETERS. (A "*" SYMBOL IS USED TO
 INDICATE THE BEST LISTED SET)

                  SET  1        SET  2        SET  3       *SET  4*
                (FEASIBLE)    (FEASIBLE)    (FEASIBLE)    (FEASIBLE)  
 ANGLE   (SV)     11.102        10.943        11.290        10.918    
 L       (DV)     300.00        295.75        305.01        295.08    
 INTENS2 (OBJ)    1.4379        1.4152        1.4167        1.4149    

                  SET  5        SET  6        SET  7        SET  8 
                (FEASIBLE)    (FEASIBLE)    (FEASIBLE)    (FEASIBLE)  
 ANGLE   (SV)     10.914        11.272        11.488        10.935    
 L       (DV)     294.97        304.53        310.31        295.55    
 INTENS2 (OBJ)    1.4150        1.4163        1.4177        1.4151    

                  SET  9        SET 10        SET 11        SET 12 
                (FEASIBLE)    (FEASIBLE)    (FEASIBLE)    (FEASIBLE)  
 ANGLE   (SV)     10.916        10.800        10.727        10.751    
 L       (DV)     295.05        291.93        290.00        290.64    
 INTENS2 (OBJ)    1.4150        1.4185        1.4181        1.4183    

                  SET 13        SET 14        SET 15        SET 16 
                (FEASIBLE)    (FEASIBLE)    (FEASIBLE)    (FEASIBLE)  
 ANGLE   (SV)     10.941        11.279        11.493        10.978    
 L       (DV)     295.70        304.71        310.42        296.70    
 INTENS2 (OBJ)    1.4152        1.4164        1.4177        1.4154    

                  SET 17        SET 18        SET 19        SET 20 
                (FEASIBLE)    (FEASIBLE)    (FEASIBLE)    (FEASIBLE)  
 ANGLE   (SV)     10.919        10.915        10.910        10.797    
 L       (DV)     295.12        295.01        294.87        291.86    
 INTENS2 (OBJ)    1.4149        1.4150        1.4150        1.4185    

                  SET 21 
                (FEASIBLE)  
 ANGLE   (SV)     10.727    
 L       (DV)     290.00    
 INTENS2 (OBJ)    1.4181